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Abstract: We design an all-dielectric Liineburg lens as an adiabatic 
space-variant lattice explicitly accounting for finite film thickness. We 
describe an all-analytical approach to compensate for the finite height of 
subwavelength dielectric structures in the pass-band regime. This method 
calculates the effective refractive index of the infinite-height lattice from 
effective medium theory, then embeds a medium of the same effective index 
into a slab waveguide of finite height and uses the waveguide dispersion 
diagram to calculate a new effective index. The results are compared with 
the conventional numerical treatment - a direct band diagram calculation, 
using a modified three-dimensional lattice with the superstrate and substrate 
included in the cell geometry. We show that the analytical results are in good 
agreement with the numerical ones, and the performance of the thin-film 
Liineburg lens is quite different than the estimates obtained assuming 
infinite height. 

© 2011 Optical Society of America 

OCIS codes: (050.6624) Subwavelength structures; (310.0310) Thin films; (230.7400) Waveg- 
uides, slab; (110.2760) Gradient-index lenses. 



References and links 

1. J. C. Maxwell, Cambridge and Dublin Math. J. 8, 188 (1854). 

2. R. K. Liineburg, Mathematical Theory of Optics (Brown U.P., Providence, 1944). 

3. J. E. Eaton, An Extension of the Luneburg-Type Lenses (Rep. No. 41 10, Naval Res. Lab., 1953). 

4. U. Leonhardt, "Optical Conformal Mapping," Science 312, 1777-1780 (2006). 

5. J. B. Pendry, D. Schurig, and D. R. Smith, "Controlling electromagnetic fields," Science 312, 1780-1782 (2006). 

6. J. Valentine, J. Li, T. Zentgraf, G. Bartal, and X. Zhang, "An optical cloak made of dielectrics," Nat. Mater. 8, 
568 (2009). 

7. D. H. Spadoti, L. H. Gabrielli, C. B. Poitras, and M. Lipson, "Focusing light in a curved-space," Opt. Express 
18,3181-3186 (2010). 

8. B. Vasic, G. Isic, R. Gajic, and K. Hingerl, "Controlling electromagnetic fields with graded photonic crystals in 
metamaterial regime," Opt. Express 18, 20321-20333 (2010). 



9. E. Hecht, Optics (4th Edition, Section 6.4) (Addison- Wesley, 2002) 

10. G. R. Schmidt, Compound Optical Arrays and Polymer Tapered Gradient Index Lenses (PhD Thesis) (University 
of Rochester, 2009) 

11. Manufacturable Gradient Index Optics (M-GRIN) (Broad Agency Announcement, Defense Advanced Research 
Projects Agency, 2010) 

12. Y. Jiao, S. Fan, and D. A. B. Miller, "Designing for beam propagation in periodic and nonperiodic photonic 
nanostructures: Extended Hamiltonian method," Phys. Rev. E. 70, 036612 (2004). 

13. P. S. J. Russel and T. A. Birks, "Hamiltonian optics of nonuniform photonic crystals," J. Lightwave Technol. 17, 
1982-1988 (1999). 

14. S. Takahashi, C. Chang, S. Y. Yang, and G. Barbastathis, "Design and fabrication of dielectric nanostructured 
Luneburg lens in optical frequencies," in Optical MEMS and Nanophotonics, (IEEE Photonics Society, 2010), 
Paper Thl-1, pp. 177-178. 

15. S. Takahashi, Design and Fabrication of micro- and nano- dielectric structures for imaging and focusing at opti- 
cal frequencies (PhD Thesis) (Massachusetts Institute of Technology, 2011) 

16. A. Vakil and N. Engheta, "Transformation Optics Using Graphene," Science 332, 1291 (2011). 

17. H. Gabrielli, J. Cardenas, C. B. Poitras, and M. Lipson, "Silicon nanostructure cloak operating at optical frequen- 
cies," Nat. Photonics 3, 461-463 (2009). 

18. T. Zentgraf, J. Valentine, N. Tapia, J. Li, and X. Zhang, "An Optical "Janus" Device for Integrated Photonics," 
Adv. Mater. 22, 2561-2564 (2010). 

19. K. K. Y. Lee, Y. Avniel, and S. G. Johnson, "Design strategies and rigorous conditions for single-polarization 
single-mode waveguides," Opt. Express 16, 15170-15184 (2008). 

20. R. Ulrich and M. Tacke, "Submillimeter waveguiding on periodic metal structure," Appl. Phys. Lett. 22, 251-253 
(1973). 

21. R. D. Meade, A. Devenyi, J. Joannopoulos, O. Alerhand, D. Smith, and K. Kash, "Novel applications of photonic 
band gap materials: Low-loss bends and high Q cavities," J. Appl. Phys. 75, 4753-4755 (1994). 

22. S. Fan, P. R. Villeneuve, J. Joannopoulos, and E. Schubert, "High extraction efficiency of spontaneous emission 
from slabs of photonic crystals," Phys. Rev. Lett. 78, 3294-3297 (1997). 

23. S. G. Johnson, S. Fan, P. R. Villeneuve, J. D. Joannopoulos, and L. A. Kolodziejski, "Guided modes in photonic 
crystal slabs," Phys. Rev. B 60, 5751 (1999). 

24. M. Qiu, "Effective index method for heterostructure-slab-waveguide-based two-dimensional photonic crystals," 
Appl. Phys. Lett. 81, 1163 (2002). 

25. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystals: Molding the flow of light, 
2nd Edition (Princeton U. P., Princeton, 2008) 

26. J. Witzens, M. Loncar, and A. Scherer, "Self-coUimation in planar photonic crystals," Selected Topics in Quantum 
Electronics, IEEE Journal of 8, 1246-1257 (2002). 

27. M. Ahmadlou, M. Kamarei, and M. H. Sheikhi, "Negative refraction and focusing analysis in a left-handed 
material slab and realization with a 3D photonic crystal structure," Journal of Optics A: Pure and Applied Optics 
8, 199 (2006). 

28. J. Zhang, L. Liu, Y. Luo, S. Zhang, and N. A. Mortensen, "Homogeneous optical cloak constructed with uniform 
layered structures," Opt. Express 19, 8625-8631 (2011). 

29. X. Chen, T. M. Grzegorczyk, B. I. Wu, J. Pacheco Jr, and J. A. Kong, "Robust method to retrieve the constitutive 
effective parameters of metamaterials," Phys. Rev. E 70, 016608 (2004). 

30. M. Hammer and O. V. Ivanova, "Effective index approximations of photonic crystal slabs: a 2-to-l-D assess- 
ment," Opt. Quant. Electron. 41, 267-283 (2009) 

31. S. M. Rytov, "Electromagnetic properties of a finely stratified medium," Sov. Phys. JETP 2, 466C475 (1956) 

32. R. Brauer and O. Bryngdahl, "Design of antireflection gratings with approximate and rigorous methods," Appl. 
Opt. 33, 7875-7882 (1994). 

33. W. Yu, T. Konishi, T. Hamamoto, H. Toyota, T. Yotsuya, and Y. Ichioka, "Polarization-multiplexed diffractive 
optical elements fabricated by subwavelength structures," Appl. Opt. 41, 96-100 (2002). 

34. J. A. Kong, Electromagnetic wave theory (EMW Publishing, 2008). 

35. S. Johnson and J. Joannopoulos, "Block-iterative frequency-domain methods for Maxwell's equations in a 
planewave basis," Opt. Express 8, 173-190 (2001). 

36. H. Gao, S. Takahashi, L. Tian, and G. Barbastathis, "Aperiodic subwavelength Luneburg lens with nonlinear Kerr 
effect compensation," Opt. Express 19, 2257-2265 (2011). 

37. A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. Joannopoulos, and S. G. Johnson, "MEEP: A flexible 
free-software package for electromagnetic simulations by the FDTD method," Comput. Phys. Commun. 181, 
687-702 (2010). 



1. Introduction 



Gradient Index (GRIN) media have been known to offer rich possibiHties for Hght manipu- 
lation since at least Maxwell's time lUl. More recent significant examples are the Liineburg 
lens ||2l, the Eaton lens O, and the plethora of imaging and cloaking configurations devised 
recently using conformal maps and transformation optics ||4]-[8l . GRIN optics are of course also 
commercially available, but the achievable refractive index profiles ^(r) are limited generally 
to parabolic in the lateral coordinates or to axial without any lateral dependence . There is 
an ongoing effort to achieve more general distributions using stacking of photo-exposed poly- 
mers ifTOllTTIl . 

For optics-on-a-chip or integrated optics applications, it is possible to emulate an effective 
index distribution n(r) by patterning a substrate with sub wavelength structures. If these are 
sufficiently smaller than the wavelength, to a good approximation they can be thought of as a 
continuum where the effective index is determined by the pattern geometry. For example, one 
can create a lattice of alternating dielectric-air with slowly varying period and fixed duty cycle, 
or with fixed period but slowly varying duty cycle (EKTSl • 

If the critical length of the variation is slow enough compared to the lattice constant that the 
adiabatic approximation is valid, the lattice dispersion diagram can be used to estimate the local 
effective index 11211131 . Refractive indices computed using a 2D approximation are valid for 2D 
adiabatically variant metamaterials where the height in the 3^^ dimension is much larger than 
the wavelength so the assumption of infinite height can be justified. According to this, we have 
designed a sub wavelength aperiodic nano structured Liineburg lens [14. 15 ] . This lens mimics a 
GRIN element with refractive index distribution n{p) = noy^2 — (p/R)^ {0 < p < R), where 
no is the ambient index outside the lens region, R is the radius of the lens region and p is the 
radial polar coordinate with the lens region as origin. The Liineburg lens focuses an incoming 
plane wave from any arbitrary direction to a geometrically perfect focal point at the opposite 
edge of the lens ll2lfT6ll. 

However, most such adiabatically variant structures are fabricated by etching holes or rods 
on a thin silicon film, whose height is less than even the optical wavelength |[6l|T4l[T5l[T7l[T8l . 
Hence, the infinite height assumption becomes questionable. Moreover, the structures are asym- 
metric since typically beneath the structure there is a substrate such as glass, whereas above the 
structure is air. Asymmetry also induces a long- wavelength cutoff in the guided modes |19]; 
therefore the thin- film metamaterial should operate in an intermediate regime where the wave- 
length is neither too large nor too small. The problem of asymmetry and finite height have been 
acknowledged in the literature on photonic crystals ll20l - [24ll . where the most common solu- 
tion is to compute a full 3D band diagram 1 25 1 . Most of them focus on photonic crystal slabs 
operating at wavelengths comparable to the periodicity, discussing phenomena such as super- 
collimation |26|, negative refraction 1271 . etc. To the best of our knowledge, the same problem 
has received insufficient attention in the context of 2D dielectric periodic or aperiodic metama- 
terial devices, especially those operating at the propagation regime of the band diagram. It has 
been briefly mentioned in |[6l[28l without giving a detailed solution. 

In our fabricated Liineburg lens design, thin-film problem is obvious where the experimen- 
tal results show dislocated and aberrated focal point (MKTS . In this paper we re-designed the 
Liineburg lens to include the finite film thickness, improving the estimate of the expected focal 
point position. To design such a lens, first we need a method for estimating effective refractive 
index of thin-film metamaterials. Several methods have been proposed in the literature. A con- 
ventional numerical approach (we refer to it as Direct Band Diagram, DBD) in photonic crystals 
derives a 3D lattice cell from the original 2D cell by surrounding a finite-height rod with large 
spaces of air above and glass substrate below ll25l . Another method takes one unit cell and re- 
trieve the refractive index by its reflection and refraction properties 1291 . These methods yield 




Fig. 1 : (a) Finite height rod lattice structure investigated in this paper, (b) 2D rod lattice structure 
assuming infinite height. 

accurate results but require either 3D band or finite-difference calculations. More heuristic (but 
faster) effective-index methods estimate a slab-waveguide effective index first and then use it 
to compute a 2D band diagram or effective index 1301 . They are generally suitable for struc- 
tures with etched substrates. In contrast, our proposal essentially reverses the order of these 
steps: we compute an effective index from the 2D cross-section first, and then incorporate it 
into a slab-waveguide mode. This is more suitable to the metamaterial regime. 

In particular, we propose the following all-analytical method for effective refractive index 
calculation. First, we replace the rods with a continuum of a certain effective permittivity £^^. 
We calculate e^^ from 2D lattice of infinite-height rods using second-order effective medium 
theory, and then substitute e^^ as the permittivity of a slab of finite thickness, acting as an 
effective guiding medium, sandwiched between semi-infinite spaces of air above and glass 
below. The geometry then becomes one of a weakly-guiding waveguide due to the small height 
of the effective guiding medium (See Fig. [2]). This weakly-guiding effect modifies the real 
part of the horizontal wave-vector component, and thus a new effective permittivity for 
the finite slab of rods is derived from the waveguide dispersion relationship. We refer to this 
method as Effective Guiding Medium (EGM). Comparing with rigorous 3D calculations, our 
method provides more physical insights, and is generally faster to compute. 

To validate our method, we compare it with the DBD method. It is shown that the results of 
both methods are in good agreement. 

2. Analytical method for effective refractive index estimation 

In this paper, without loss of generality, we investigate a silica glass slab covered by a square 
lattice (lattice constant a = 258 nm) of silicon rods of finite height h = 320 nm, variable radius 
r (0 < r < a/V2) and immersed in air, as illustrated in Fig. [Ha). The free space wavelength 
of light is chosen as A = 6a = 1550 nm. This choice of a is small enough to insure that we 
remain in the metamaterial regime and in the propagating regime of the band diagram; and large 
enough that the rods can be accurately fabricated by nano-lithography |[l4l[T5l and we do not 
reach the long-wavelength cutoff regime for the asymmetric waveguide, as mentioned above. 
The dielectric permittivity constants for glass and silicon are Cgiass =2.25 and ^silicon = 12.0, 
respectively. These media are non-magnetic, so the relative permeability is taken as /i = 1 
throughout this paper. The glass slab height is assumed to be much larger than the height of the 
rods and the free space wavelength of the light. The corresponding 2D structure with infinite 
height rods and without glass substrate is shown in Fig. [Hb). We now proceed to describe 
all-analytical method, EGM, for analyzing these two geometries. 



2.1. Effective guiding medium (EGM) method 

The EGM method requires analysis of a three-layer structure: (I) air, (II) effective medium 
waveguide and (III) glass, as shown in Fig. [21 The effective permittivity of the guiding medium 
is calculated from the second-order effective medium theory in 2D which have been derived by 
various authors ||3T][33. This theory starts from the effective refractive index of ID sub wave- 
length grating composed of air and dielectric with index n. Under TE (electric field parallel to 
the grooves) and TM (electric field vertical to the grooves) polarization incidence the effective 
index can be summarized, respectively, as 



4e = «oTE+Y(^)/(l-/)'(«'-l)^ (1) 

4m = «otm + y(^)/(1-/)'«otm«ote(;|2-1 ) ' (2) 



«OTE =/«' + (!-/), «0TM = l/(4 + (l-/)) (3) 



where 



are the zeroth-order effective refractive indices, T is the period of the grating and / is the 
filling factor of the dielectric grooves. The effective indices of corresponding 2D sub wavelength 
structures are then estimated as a combination of ID structures I 



^2D-TE = y^-f^fK^^ (4) 



«20-TM = (Va-/)+MM+W ,.^(,:^^)^J /2 (5) 

for both TE and TM polarizations. Note that TE and TM polarizations mentioned in this paper 
are an approximation since the fields are not purely polarized in 3D structures. A more exact 
way to describe them is TE-like/TM-like, where electrical field is mostly parallel/vertical to the 
grooves ||25]| . However, this is still an approximation because the waveguide is asymmetric so 
there is no horizontal mirror symmetric plane. The second-order terms used in Eq.[TJ&|2] better 
approximate the effective index in the case that the wavelength is not very large comparing with 
size of unit cell, e.g. X = 6a used in this paper. Most current metamaterial device designs are 
using the zeroth-order approximation only O, even when the unit cell size is not far smaller 
than the operational wavelength. This is fine for those devices where high accuracy results 
are not important. However, for devices such as Liineburg lens, all waves are focusing to a 
single point so light manipulation is more challenging. Therefore, more precise effective index 
prediction is needed and second-order corrections are included. 

The dispersion relation of the effective guiding medium, i.e. the relationship between and 
(0, is governed by the guidance condition of an asymmetric dielectric waveguide for both TE 
and TM polarizations ll34l 



(TE :) i^n(k^yh) = ji 2 rr? 9/9 //? 9 = ^TEihyh), (6) 



(TM :) tan(^ii^/z) = 



Hy 



Fjuihyh), (7) 



Fig. 2: Effective guiding medium (EGM) approximation of 2D finite height rod lattice structure. 
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Fig. 3: Graphical solutions of wave guidance condition (Eq. [6]-[7]) for TE (a) and TM (b) po- 
larizations. Blue and red lines are the left and right hand sides of these equations, respec- 
tively. Operating frequencies = 0.11 x Inc/a, 0)2 = 0.16 x Incja, = 0.14 x Inc/a and 
0)4 = 0.18 X iKc/a. Rod radius r = 0.50a. 



where = £j[CO^ / — k^^ is the phase-matched propagation constant. These equations can 
be solved by a graphical method and an example is illustrated in Fig.O It is observed that one 
and only one intersection is obtained for each frequency, meaning that only one fundamental 
mode is supported. Full dispersion relations kz{co) are shown in the following section. 

The EGM method described above is compared with the conventional DBD method. To 
apply the DBD method, we need to calculate the band diagram of the 3D super cell shown in 
Fig-ffla). The supercell height is taken as large sls H = 20a to better emulate the real structure 
of Fig. [Ha), where the air and glass spaces tend to infinity. In other words, we seek to minimize 
the interference between neighboring unit cells along the vertical (y) direction. We used the 
MIT Photonic-Bands (MPB) mode solver ||35]| to calculate the dispersion diagram. In Fig.|4tb- 
c) we show an example MPB result for our chosen lattice and the specific value r = 0.5a, for 
temporal frequency ft) = 1/6 x Inc/a. From Fig.^b) we observe that for the chosen values of r 
and ft), the isofrequency contour 1251 is almost a circle, indicating that this unit cell is isotropic. 
Therefore, when using DBD in this particular geometry, it is sufficient to consider kz{(o) only. 
However, this is not generally true in other geometries as r or ft) increase. 

Figure Hfc) shows the mode shape for the same geometry. It can be seen that the field is 
effectively concentrated near the silicon rod portion of the cell. The relative intensities at two 
horizontal cell boundaries y = ±///2 were 5.6 x 10~^ and 3.8 x 10~^ at the top and bottom, 
respectively, compared to the peak value that occurred at j = 159 nm from the rod base. This 



validates our choice of H as sufficiently large. 




(a) (b) (c) 

Fig. 4: (a) The supercell used in the DBD method for the finite height rod lattice structure, 
(b) Isofrequency contour of the supercell with r = 0.50a where the first band only is shown. 
Labels on the lines denote the corresponding normalized frequency (Oa/lnc. The bold blue line 
corresponds to the wavelength X=6a used in this paper, (c) Field distribution of the waveguide 
slab at a particular x slice. Color shading denotes magnetic field {Hy) distribution and black 
contours illustrate silicon rods. 

Comparing with the DBD method, the EGM method can provide deeper physical insights 
with all-analytical solutions, and is generally faster since it avoids solving numerical electro- 
magnetic solutions in 3D. 

2.2. Effective refractive index and rod radius relationship 

In this section, the relationship between the effective refractive index and rod radius is calcu- 
lated. The results of EGM method are compared with the ones obtained from DBD method. 

Figure [5ta) shows the dispersion relation of the finite-height rod lattice calculated with both 
DBD and EGM methods, as well as with the 2D (infinite rod height) assumption, for rod radius 
r = 0.5a. Based on the dispersion relation, effective refractive indices for unit cells with differ- 
ent rod radii are calculated as ^gff = ck^/ co, shown in Fig. [3b). The results given by the DBD 
and EGM methods are in good agreement with each other, with maximum percentage errors 
of 7.3% and 6.0% for 2D and 3D cases, respectively. It is observed that the effective refractive 
indices of the finite-height rods are significantly different than those assuming infinite height. 
This is to be expected due to weak guidance: as can been seen in Fig. He), a large portion of 
the field extends outside the rods to spaces of air and substrate. When the rod radii are below 
certain values (0.17a for TE and 0.35a for TM), the propagation modes are not guided so the 
effective indices are not shown. The discontinuities observed in the 2D effective index curves 
for DBD method beyond certain values of rod radii (0.40a for TE and 0.49a for TM) result 
from the emergence of a photonic crystal bandgap at these values. At this frequency range, 
even though the 2D infinite-height lattice is within the bandgap, the confined (slab waveguide) 
geometry is still propagating; this is because the light is mostly outside the dielectric region, so 
propagation takes place in the free space (hence the lower index). To calculate the propagation 
constant in this regime, we still need an effective index value and EGM provides it (it turns out 
to be large than 3, typically). 
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Fig. 5: (a) Comparison between the dispersion relation for finite-height siHcon rod lattice 
[Fig-ffia)] calculated from the EGM and DBD method, and the dispersion relation for infinite- 
height 2D rod lattice [Fig. [Tib)]. For each case, the two lowest bands representing the TM and 
TE modes are shown, (b) Relationship between effective refractive index and rod radius cal- 
culated from both methods, compared with the relationship for infinite-height 2D rod lattice. 
Free space wavelength of light is A = 6a = 1550 nm. 

3. Optimal design of the subwavelength Luneburg lens 

We re-design and numerically verify the subwavelength Liineburg lens 1121141151361 . which was 
previously designed under 2D assumption. Here, we still design the Luneburg lens as a structure 
consisting of finite-height rods with adiabatically changing radius r across the lattice of fixed 
constant a. At each coordinate p, we emulate the Luneburg distribution n{p) = ^o\/2 — {p/R)^ 
by choosing the rod radius r at coordinate p from Fig. Ob) such that = ^(p), as opposed to 
using = n{p). The design has to be carried out separately for the TE and TM polarizations. 
The ambient index is chosen as = 1-53. 

Figure illustrates the lens structures and the corresponding 3D finite-difference time- 
domain (FDTD) simulation results for the actual adiabatically variant thin-film nano structured 
Luneburg lens performed by MIT Electromagnetic Equation Propagation (MEEP) 1371 . With 
plane wave illumination, almost diffraction-limited focal points at the edge can be observed for 
both TE and TM polarizations. For a more computationally efficient and intuitive representa- 
tion we also ray-traced the field inside the Luneburg structure using the adiabatic Hamiltonian 
method |[l2l[T3l . The ray position q and momentum p are obtained by solving the two sets of 
coupled ordinary differential equations 



where if (q,p) = ()L)(p,k) is obtained from the dispersion diagram at each coordinate |q| = p 
and for k = p. Ray tracing results are superimposed in Fig. [6l with FDTD results, and are seen 
to be in good agreement. Furthermore, as a comparison, similar thin-film Luneburg lens is 
designed using the DBD method and simulation results are shown in Fig. [71 It is observed that 
results of the all-analytical EGM method design agree with those from the DBD method. 

To compare the redesigned lens (3D, finite height) with the original design (2D, infinite 
height), we repeated the design using the values of refractive indices predicted by the dis- 
persion relation of the infinite-height rod lattice (see Fig. Ob) blue and red solid curves). In 
this case, we are forced to use TM polarization only because the TE polarization reaches the 
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Fig. 6: (a) Top view and side view of the thin-film sub wavelength Liineburg lens designed by 
EGM method for TE mode and (b) the corresponding 3D FDTD and Hamiltonian ray tracing 
results, (c) Top view and side view for TM mode and (d) the corresponding 3D FDTD and 
ray tracing results. Red circles outline the edge of Liineburg lens, where radius 7? = 30a. Blue 
lines are the ray tracing results and color shading denotes the field [^3; for (b) and Ey for (d)] 
distribution, where red is positive and blue is negative. 
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(c) (d) 

Fig. 7: Structure and the corresponding 3D FDTD and Hamiltonian ray tracing for the thin-film 
subwavelength Liineburg lens shown in Fig.[6l but designed by the DBD method instead. 

bandgap for relatively small value of r, not leaving enough room to implement the Liineburg 
profile with rod radius r large enough to be robust to practical lithography and etching methods 
(in our experiment, this requires r > 0.27a II14II15I ). It can be observed from the FDTD and 
Hamiltonian ray-tracing results shown in Fig. [S] that the focal point is outside the lens edge and 
it is strongly aberrated. This is in good agreement with the experimental results of the original 
design Cma. 
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Fig. 8: FDTD and Hamiltonian ray-tracing results of the sub wavelength Liineburg lens made 
of finite height silicon rods, but designed assuming infinite height. The color conventions are 
the same as in Fig.Ob&d) 



